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Abstract 

Prof. Raman found that the Indian drums ah produced harmonics, 
whereas the western drums do not, a fact that is proved theoreticaUy. 
He also found that that the density distribution present as a central 
black patch on the Indian drums is not of the type r~". However, 
R. Siddharthan, P.Chatterjee and Vikram Tripathi have stated in their 
paper |^ that they found the harmonics of the Indian Drums all in order, 
except the fundamental, which was higher than what it should have 
been. This article presents two theoretical models (or distributions) one 
of which was found to yield frequency ratios identical to those found by 
R. Siddharthan et al. 



^see Physics Education Oct. 1994. 




Figure 1: A view of the Tabla (photo courtesy of TABLA.COM, Copyright 
1999). This article is concerned with the right-hand side tabla, which can be 
seen in this figure on the left. 
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1 Introduction - Music, Harmonics and Over- 
tones 



Sound is caused by alternate compressions and rarefactions of pressure in the 
air. The ear-drum vibrates with these so-called longitudinal waves, and that 
is how the sense of sound is had. Any arbitrary sound, however, does not 
please the ear. The normal human ear can, in general, differentiate between 
what is pleasing and what is noise. This section gives the important differences 
between Music and Noise. 

1.1 Music 

Music has its own definite wave patterns, so that any arbitrary sound cannot 
be called music. One can, by looking at the wave patterns, find out whether 
a given sound is musical or not. Music has the following two attributes that 
distinguish it from noise: 

1. Pitch: The pitch of a wave pattern characterises the time after which 
the wave pattern repeats itself. For the sensation of pitch, therefore, it is 
essential that the waveform be periodic. A non-periodic waveform 
will not give the sensation of pitch, and can hence be labeled noise. 

2. Timbre: Timbre refers to the quality of sound. A tuning fork, a sitar, 
a tanpura and a piano may all have the same pitch, yet their timbre may 
be completely different. Pitch refers to the frequency of repetition of 
the waveform, whereas timbre depends on the shape of the waveform. 
What timbre really refers to is the relative amplitudes of the various 
component single frequencies that make up a particular sound wave. 

1.2 Overtones and Harmonics 

Any complex waveform may be expressed as a linear superposition of a number 
of sine waves of different frequencies. Notes produced by musical instruments 
consist of sine waves of one fundamental frequency and of higher frequencies 
called overtones. For the resultant waveform to be periodic, the overtones 
have to be integral multiples of the fundamental frequency. Overtones that 
are integral multiples of the fundamental frequency are called harmonics. 

Unless a majority of overtones in a particular sound are harmonic, the 
waveform will not be periodic and will not have a discernible pitch. Thus 
for a note to sound musical 0, a majority of the overtones must be 
harmonic. 

"^Musical., as used here just means having a definite pitch, which is very important in 
Indian Classical accompaniments; however, this does not mean that the Western drums are 
'non-musical', they produce rhythm, just like their Indian counterparts, see section 1.3. 
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1.3 Indian and Western Drums 



The Tabla is the most well-known of all Indian drums. Amir Khusro created it 
when he split the ancient Indian drum, the Pakhawaj into two parts. Of these, 
the right drum has a black patch in the centre, as mentioned earlier. This is 
made of a mixture of iron, iron oxides, resin, gum etc. and is stuck firmly on 
to the membrane. The thickness of this patch decreases radially outwards. 

The left hand drum has a wider membrane, and has a black patch similar 
to the one in the right drum, except that it is unsymmetrically placed on one 
side of the membrane. It is worth noting here that the left hand drum is 
not used to produce harmonics but to provide lower frequencies in the overall 
sound while the Tabla is being played. 

The key difference between Indian and western drums is the absence of the 
central loading in the case of the Western drums. It is shown mathematically 
in section 4 that a uniform circular membrane cannot produce harmonics. 
In section 5, we then investigate how the drum may be made harmonic by 
considering two theoretical models, i.e. radial density distributions. 

Thus, this is a study of how the density variation of the membrane 
affects the frequencies of the overtones. 

Indian Classical music is such that a recital consists (usually) of a single 
performer, and a couple of accompaniments. In this situation, the accom- 
paniment (usually a drum, the Tabla for instance, must be harmonic, since 
otherwise, the aharmonicity would completely disrupt the recital. 

Western drums, however, play a very different part in Western Classical 
music. They provide a rhythm to the music produced by the rest of the 
instruments, and (strictly speaking) do not need to produce harmonics. 



2 The Wave Equation and its solution 

The wave equation, in the case of a membrane, describes the relationship 
between the displacement of a point on the membrane changes with position 
and with time, when it is disturbed in some way. 
The general form of the Wave equation is : 

For a circular membrane, in polar co-ordinates, the wave equation may be 
written as 

d'^u Idu 1 d'^u 1 d'^u 

where (r, 0) describes any arbitrary point on the membrane that is assumed to 
have: 

1. Uniform mass per unit area 

2. Uniform tension per unit length 
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and u{r, 6, t) is the vertical displacement of the point at a time t. Also, 

- ^ (3) 

where c is the speed at which the wave travels on the membrane, r is the 
tension per unit length in the membrane and a is the mass per unit area of 
the membrane. 

This equation is solved using the technique of separation of variables, where 
it is assumed that 

u{r,e,<j)) = R{r)e{9)T{t). (4) 

Here, R depends only on r, © only on 9, and T only on t. Making this 
substitution, and solving, the three quantities T, © and R are found to be: 

1. The T solution: This is a simple sine (or cosine) function, indicating 
that the membrane undergoes simple harmonic oscillations in time. Its 
general form is: 

T = cos {ut + (j)) , (5) 

where omega = ck and is a phase which depends on the initial condi- 
tions. However, it can be set to zero without loss of generality. 

2. The solution: This tells us how u changes with angle, i.e. how 
displacement changes as we move along a circular path of fixed radius. 
This part is: 

e^cos{m0 + ip). (6) 

ip is again a phase which depends on initial conditions, and may be put 
to zero without loss of generality. There are certain values of 9 for which 
0(^) reduces to zero. Thus, u would, at all times be equal to zero on all 
diameters along which m9 is an integral multiple of |. 

Nodal Diameters of a drum are those diameters that remain stationary 
while the rest of the drum is vibrating. For a particular value of m, there 
are m symmetrically placed Nodal Diameters. 

3. The R solution: This describes how the displacement of the membrane 
above the plane of rest (i.e. u) changes as we move outward, along a 
radius. There are certain circles that remain stationary on the membrane 
while the rest of the drum vibrates. This happens whenever the function 
R{r) crosses the value zero. These circles are called Nodal Circles and 
their radii correspond to those values of r at which R{r) reduces to zero. 
The R solution is given by: 

R{r) = J^ikr), (7) 
where Jm is the Bessel function of the first kind of order m. 
A few Bessel functions are plotted in Fig. 2 and Fig. 3. 
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Bessel functions of order and 1 
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Figure 2: Bessel Functions of order and 1. 

Plot of J(x) 
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Figure 3: Bessel Functions of order 3,4,5 and 6. 
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Figure 4: The first nine Normal Modes of a circular membrane. Quantities in 
brackets are {Nodal Diameters, Nodal Circles). Arrows indicate the direction 
in which the various parts, separated by Nodal circles and diameters, vibrate. 

2.1 Normal Modes 

The complete solution may be written as: 

u{r, 9,t) = Jm{kmnr)cos{m9)cos{umnt) (8) 

The most general solution is thus a superposition of various different modes, 
where in each mode, the whole drum vibrates with one frequency, has m Nodal 
Diameters and n— 1 nodal circles. These modes are called the Normal Modes 
of the drum. The first nine Normal Modes are shown in Fig. 4. 

2.2 The Boundary condition 

The membrane is bound in the form a circle over the head of the wooden 
shell. Therefore, the periphery is always stationary. This has to be taken into 
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account when the solution is being calculated. The condition may be stated 
mathematically as follows: 

u{r, e,t)^0 (9) 

at r = a or 

Jnikmna) = 0. (10) 

That is, the boundary of the membrane should correspond to one of the zeroes 
of the Bessel function of any order n. This will yield the values of frequencies 
that are 'allowed' by the boundary condition. However, we are interested only 
in the relative ratios of these frequencies. 

Let us quantify the 'allowed' frequencies. Let be the mth zero of the 
Bessel function of order n. Then the above boundary condition becomes: 

a 

i.e. the allowed frequencies are proportional to the zeroes of Bessel functions. 
As noted earlier, however, the values of these zeroes have no integral ratio to 
each other, so that the ratios of the frequencies are also non-integral. 

Thus, the simple circular membrane, and hence also Western 
drums, cannot produce harmonics. 



3 Making the Drum Harmonic 

Indian drums, in general, play a very different role as accompaniments in 
Classical music as compared to their western counterparts, as mentioned in 
sections 1.2 and 1.3. Since a uniform membrane does not give us harmonics, we 
tried solutions to a membrane with a density variation. The simplest possibility 
is a loading which varies only with r. We looked at the right hand tabla 
and tried various symmetric density distributions that resembled the actual 
loading. We describe here the one that was found to be successful. But first, 
we note the change produced in the wave equation. 

Recall that = ^. Now, however, p = p{r), so that c^(r) = This 
causes the wave equation to change from 

IdR^ /o^ _ _^ (12) 

for a uniform membrane to 

(PR IdR f 




+ --7-+ 37-t--:tM? = (13) 



or, equivalently. 



dr^ r dr \ c^(r) 



d^R IdR fp(r)io^ „ , 



dr'^ r dr \ T 
for a loaded membrane 
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Plot of the continuous loading 
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Figure 5: A side view of the loading, as is seen on the Tabla, represented as 
the function alog{ro — r) + b. 



The Radial equation can be solved for various distributions. We tried two 
kinds of loading patterns: 

1. The Step function: Concentric rings with varying density were con- 
sidered. This gave good results, as can be seen in the table. 

2. The Continuous Loading: The actual loading on the tabla is a step 
function to begin with. However, even though it is stuck in parts, the 
loading becomes more or less continuous after the tabla has been played 
for some time. The loading therefore looks quite like the function plotted 
in Fig. 5. 

It was also found that if an exponential function is included toward the pe- 
riphery, the results obtained for both the distributions improve. The form of 
the function used, along with the continuous loading, was ce*, where c, d are 
constants. The actual function used as loading is shown in Fig. 6. 
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r = 0.62 r = 0.93 r=1 



Figure 6: The actual function used. 

4 Calculation and results 

Both the distributions given in the previous section were put into the Radial 
part of the wave equation, and the equation was solved numerically using the 
second-order Runge-Kutta method. 

In writing down the wave equation for the loaded membrane, the following 
two assumptions were involved: 

1. Normal Modes exist even in the loaded membrane 

2. Tension per unit length is the same throughout the membrane, even in 
the loaded part. 

4.1 Calculation 

It may be seen from figure that the function chosen for Continuous loading 
varies slowly with r near the centre of the membrane. Also, in the discrete 
loading, the density remains constant in a particular concentric circle, so that 
that the density is not varying in the centre. Due to this reason, the initial 



10 



conditions (for starting the solution of the wave equation by the Runge-Kutta 
method) were assumed to be identical to those of the Bessel functions. 

The allowed frequencies were found in the following way: First, the Radial 
equation was written as 



This is the equivalent of the Bessel equation of order m for the loaded case. 

Now, keeping m = 0, the value of k' was varied and the solution plotted as 
a graph on the screen until its value became zero at the boundary. The corre- 
sponding value of k' was noted, k' was then increased, until the solution again 
became zero, but this time the solution passed through zero once, meaning 
that there was a nodal circle. The process was continued for a particular value 
of m. Then, m was increased by one, and the same process was continued, 
and all the allowed frequency values were noted. 

4.2 Problems with the calculation 

The method mentioned in the previous section works fine for m = and rn = 1, 
since the initial values of the Bessel functions and their derivatives are non-zero 
for the zeroth and first orders. However, for m = 2 and higher orders, values 
of both the Bessel function and its derivative become zero at r = 0. This leads 
to the solution becoming zero at all points for m = 2 and higher orders. 

This happens because the Runge-Kutta method depends on the initial value 
of the solution and its derivative (i.e. the value at r = 0). Since the iteration 
begins with the derivative and the initial value as zero, it continues to be so 
for further values of r. 

To get around this problem, iteration was begun not from r = but from 
an infinitesimally small value, in this case r = 0.0001. The values of the 
Bessel function and its derivative at this point were computed from the series 
expansion of the Bessel function, considering the first four terms. The rest of 
the procedure was as in the previous section, but this 'initial value' had to be 
calculated for each order separately. 

4.3 Results 

The most surprising result we got was that for the continuous loading, all the 
overtones were found to be nearly harmonic, but the "fundamental" itself was 
higher than what it should have been. In particular, the ratios were found to 
be 1.07:2:2:3 etc. These results are summarized in Table l.The base in column 
3 is the 2nd normal mode, with 1 nodal diameter. This is done to illustrate 




(15) 



where 




(16) 



11 



TTn n 1 A/f nfl p ^ 

l\Ul 1 1 uUjI IvIUU/V^ 


Frequency ratio 




Unloaded 


Continuous Loading 


Multiple rings 


Tabla 


1 (0 0) 


l.UU 


i.O/ 


i.OO 


1 nn 
1.00 




i.oy 


z.UU 


i.y6 


o 
z 




0.i4 


o no 
z.yo 


O CiQ 

z.yo 


Q 
O 




o on 


o on 


Q no 


Q 
O 




o.oo 


4. (JO 


/I no 


A 

4 




o no 


4.UU 


o.yo 


4 




O 1 


O.Oi 


o.Oz 


r 





o.oU 


o.Ui 


o.UU 


r 





3.60 


5.02 


4.80 


5 


10(1,2) 


4.24 


6.02 


5.20 




11(1,3) 


5.55 


7.80 


7.03 




12(2,2) 


4.85 


7.00 


5.90 




13(3,1) 


4.06 


6.04 


6.02 




14(4,1) 


4.60 


7.09 


7.05 





the fact that the fundamental is absent and the lowest eigenvalue is 1.07 times 
the fundamental. 

Thus, a theoretical model with the first seven harmonics is possible, which 
also resembles the actual loading pattern on the Indian drums. However, with 
the kind of loading that the Indian drums have, the fundamental is absent 
and a slightly higher note is present instead. This is in accordance with recent 
observations of the frequencies present in the Indian drums. 

5 Discussion 

The results table shows certain modes with relative frequency ratios that differ 
from the exact harmonic ratios by small amounts (e.g. 0.01, 0.02). However, 
the minimum difference required between two frequencies so that a normal 
human ear can distinguish between them is 6 — IHz. Considering that the 
fundamental note in Indian Classical {Sa of the base Saptak) is 24:0Hz, these 
differences work out to 2AHz to 5Hz. Thus, the human ear is not able to 
make out these differences. 
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